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Money in Gas-Like Markets: Gibbs and Pareto Laws 
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We consider the ideal-gas models of trading markets, where each agent is identified with a gas 
molecule and each trading as an elastic or money-conserving (two-body) collision. Unlike in the 
ideal gas, we introduce saving propensity A of agents, such that each agent saves a fraction A of its 
money and trades with the rest. We show the steady-state money or wealth distribution in a market 
is Gibbs-like for A = 0, has got a non-vanishing most-probable value for A / and Pareto-like when 
A is widely distributed among the agents. We compare these results with observations on wealth 
distributions of various countries. 
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I. INTRODUCTION 



In an earlier conference in Kolkata in 1994, many lead- 
ing Indian economists from the Indian Statistical Insti- 
tute and physicists met and discussed about the possible 
formulations of some economic problems and their solu- 
tions using tricks from physics [lj . In one of these papers 
0, possibly the first published joint paper with both 
physicist and economist Indian coauthors, the possibility 
of a kinetic theory of (ideal) gas-like model of trading in 
the market was discussed. Among other things, it tried to 
identify, from the known effects of various fiscal policies, 
equivalence of the (kinetic) energy of the gas molecules 
(money) and the temperature (average money in the mar- 
ket). Such a "finite temperature" market model and 
the corresponding distributions were also noted by oth- 
ers 0, H| • With the possibility of putting more than one 
agent in the same (micro) state, identified by the price or 
money income of the agent in the market, the likely dis- 
tribution was concluded there to be Bose-Einstein like, 
rather than Gibbs like . These studies of course had the 
limitation of absence of any comparison with real income 
distributions (in any market or country) . In a recent pa- 
per by Dragulescu and Yakovenko a simple (trading) 
market model was developed with fixed (total) money 
and number of agents in the market. Random two-agent 
exchanges (with local money conservation) lead to Gibbs- 
like steady income distribution. This was also confirmed 
by simple numerical simulations. Modifications due to 
savings was studied simultaneously @. In a very recent 
review Q , a popular introduction to these developments 
is given. 

Saving propensity among the agents, a very selfish and 
local feature of the tradings, introduce in effect some 
global co-operative feature (cf. We show that a 

fixed and uniform saving propensity of all the agents in 



the market shifts the most-probable money of the distri- 
bution away from zero (as given by Gibbs for zero sav- 
ings), while a random distribution of saving propensity 
among the population can give the Pareto (power) law 




P(m) 



-(!+„) 



(1) 



for the wealth or money (to) distribution. We intend to 
discuss here in brief the effects of various kinds of savings 
on the ideal gas-like money distributions in the above- 
mentioned market models, and compare our observations 
with those from real markets. 



II. AN IDEAL GAS-LIKE MARKET MODEL 

Let us consider a simple model of a closed economic 
system where the total amount of money M and the to- 
tal number N of agents are fixed. No development (pro- 
duction) or migration (death/birth of agents) occurs and 
the only economic activity is confined to trading. Each 
agent i, individual or a corporate, possess a money m,(f) 
at (discretised) time t. Time changes after each trad- 
ing. In any trading, two randomly chosen agents % and j 
exchange their money such that their total money is (lo- 
cally) conserved and none ends up with negative money 
(debt not allowed): 



rrii(i) + TTij (t) = m. L (t + 1) + m,j(t + 1) 



(2) 



where m,(f) > for all i and t; 53;=i m i = ^d. Since 
money is conserved, in the steady state (t — » oo), the 
probability P(m) of the density of people with money m 
will satisfy 
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P(m 1 )P(m, 2 ) = P{mi + m 2 ) 



(3) 



2 



which corresponds to the Gibbs distribution 



P(m) = (1/T)exp (-m/T) ; T = M/N. (4) 

Numerical simulations in the model also confirmed the 
steady state distribution of money, no matter what initial 
distribution of money the agents had, to be Gibbs' one: 
after sufficient number of tradings, most of the agents 
end up with very little money! This result is quite ro- 
bust (and realistic too!). In fact, several variations of 
the trading, and of the 'lattice' (on which the agents can 
be put and each agent trade with its 'lattice neighbours' 
only), whether compact, fractal or small- world |3, leaves 
the distribution unchanged. Some other variations like 
random sharing of an amount 2m-2 only (not of mi + TO2) 
when mi > (trading at the level of lower economic 
class), lead to even drastic situation: all the money in 
the market drifts to one agent and the rest of the agents 
all become truely pauper [y, Q ! Attempts have also been 
made to get Pareto-like power-law distribution here 
with changed definition of entropy or the conservation 
law (cf. eqn. |J3}). 

Chakrabarti and Marjit Q argued for the Bose- 
Einstein like distribution (rather than Gibbs) in such 
a market (with the temperature T similarly identified 
with the average money M/N per agent) , as one can put 
more than one agent in the same economic state (spec- 
ified by the income) and the maximisation of the con- 
sequent entropy. For the possibility of adding and sub- 
tracting agents into/from the market, one similarly needs 
(negative) "chemical potential" which becomes zero at a 
finite temperature or money level in the market, when 
the "Bose condensed" fraction of the agents will fall out 
of the market distribution and might be identified as un- 
employed. 

The real income distributions did not indicate so far 
anything like the Bose distribution; rather considerable 
evidences support the possibility of Gibbs like distribu- 
tion J3J in the income (almost for 90% of the low-income 
range) of various countries (see e.g., [I(J], see also data in 

El; Fig. 2). 



III. MODEL WITH FIXED SAVING 
PROPENSITY OF THE AGENTS 

Here we assume |5( that each economic agent i saves 
a fraction A of its money mi(t) before the trading at 
time t. We again assume that an agent's money is non- 
negative and no debt is permitted. Let us now consider 
an arbitrary pair of agents i and j , who get engaged in a 
trade, and their money TOj(i) and rrij(t) before the trade 
change respectively to 

rrii(t + 1) = rrii(t) + Am; mj(t+l) = m,j{t) — Am (5) 



where Am = (1 — X)[e{mi(t) + mj(t)} — mj(t)]; (6) 



with e as any random fraction. As may be checked by 
straight-forward substitution, this kind of trading again 
satisfies eqn. (2), while each agent saves a fixed fraction 
A of its money before the trade and exchanges randomly 
(with fraction e) the rest of the money. 

One finds here that at A = the market becomes non- 
interacting and the steady state money distribution be- 
comes the Gibbs' one. For any nonvanishing A, the equi- 
librium distribution becomes asymmetric Gibbs-like (see 
inset of Fig. 1) with the most-probable money m p per 
agent (corresponding to the peak in P(m)) shifting away 
from m = (for A = 0) to M/N as A -> 1 @. This 
self-organising feature of the market, induced by sheer 
self-interest of saving by each agent without any global 
perspective, is very significant as the fraction of paupers 
decreases with saving fraction A and most people end up 
with the average money in the market (socialists' dream 
achieved with just people's self-interest of saving)! In- 
terestingly, self-organisation also occur in such market 
models when there is restriction in the commodity market 
|l2| . Although this fixed saving propensity does not give 
yet the Pareto-like power-law distribution, the Marko- 
vian nature of the scattering or trading processes (eqn. 
©) is lost and the system becomes co-operative. Indi- 
rectly through A, the agents get to know (start interact- 
ing with) each other and the system co-operatively self- 
organise towards a most probable distribution (m p ^ 0) . 
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FIG.l: Money distribution P(m) in the model for distributed A 
(0 < A < 1). Inset shows P(m) for three typical values in the fixed 
A case (including A = 0; Gibbs). For both cases, N = 100 play 
with average money per agent M/N = 1. 



IV. MODEL WITH RANDOM SAVING 
PROPENSITY OF THE AGENTS 

We now consider a market again with fixed N and 
M but with random saving propensity Ai (0 < Xi < 1) 
fixed or "quenched" for each agent (Xi are independent of 
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trading or t, but vary randomly from agent to agent) |13| . 
One again follows the same trading rules as mentioned 
in the previous section (eqn. QJ), except that 



Am = e(l - \j)mi(t) - (1 - Ai)(l - e)m^) (7) 

here; A, and Xj are the saving propensities of agents i and 
j. We first take a market with N agents, each having a 
fixed saving propensity A distributed independently ran- 
domly and uniformly (white) within an interval to 1. 
Having assigned each agent i the saving propensities Aj, 
and starting with an arbitrary initial (uniform or ran- 
dom) distribution of money among the agents, we start 
the tradings. At each time, two agents are randomly 
selected and the money exchange among them occurs, 
following the above mentioned scheme. We check for the 
steady state, by looking at the stability of the money 
distribution P(m) in successive Monte Carlo steps t. 

In Fig. 1, we show the money distribution P(m) vs. 
to (in units of M/N) for N = 100, M/N = 1, after 
averaging over 10 6 initial configurations (A, distribution 
among the agents) at t/N =10,000. There is an initial 
growth of P(m) from m = 0, which quickly saturates 
and then a long range of power-law decay in P(m) for 
large m values (for less than 10% of the population N 
in the market) is observed (for more than two decades 
in m). This decay, when fitted to Pareto law gives 
v = 1.03 ±0.03. 

Income tn (in thousand dollars) 




FIG. 2: Cum ulative distribution Q(m ) = P(m)dm of wealth m 
in USA 10] in 1997 and Japan 11] in 2000. Low-income group 
follow Gibbs law (shaded region) and the rest (about 5%) of the 
rich population follow Pareto law. The inset shows the cumulative 
distribution for a model market where p = 0.9 fraction of the agents 
trade randomly without any saving and the rest 1 — p fraction trades 
with their saving propensities distributed uniformly between and 
1. The dotted line (for large m values) corresponds to v = 1.0. 

We now investigate on the range of distribution of the 
saving propensities within the population. If a certain 
fraction p of the population trades in the market ran- 
domly without any saving while the rest (1 — p) fraction 
have a quenched distribution of their saving propensi- 
ties (0 < A; < 1), we observe that for larger values of p 
(p > 0.8), the distribution is Gibbs-like for low-income 
group and has a power law tail for the high-income group. 
The range of validity of Gibbs law increases with increas- 
ing p. However, the exponent v does not change with 
p. These model wealth distributions P(m) compare very 
well with the wealth distributions of various countries: 
Data suggests Gibbs like distribution in the low-income 
range [X(| (more than 90% of the population) and Pareto- 
like in the high- income range [ll| (less than 10% of the 
population) of various countries (Fig. 2). 

We also considered annealed randomness an the saving 
propensity A: here Aj for any agent i changes from one 
value to another within the range < Aj < 1, after each 
trading. Numerical studies for this annealed model did 
not show any power law behavior for P(m); rather it 
again becomes exponentially decaying on both sides of a 
most-probable value m p (A), similar to fixed A case. 
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SUMMARY AND CONCLUSIONS 



We have considered ideal-gas models of trading mar- 
kets. In these models, we introduce saving propensity A 
of agents, such that each agent saves a fraction A of its 
money and trades with the rest. We show the steady- 
state money or wealth distribution P(m) in the market 
is that of Gibbs J3J for A = 0, has got a non- vanishing 
most-probable value for A > (but fixed for all agents), 
and one gets Pareto distribution (|TJ with v ~ 1.0 when 
A is widely distributed among the agents. These results 
in simple ideal-gas like market models also compare well 
with real market observations. 
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